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We propose an algorithm that transforms a digitized color image so as to simulate for normal. observers, the 
appearance of the image for people who have dichromatic forms of color blindness. The dichromat's color con- 
fusions are deduced from colorimetry, and the residual hues in the transformed image are derived from the 
reports of unilateral dichromats described in the literature. We represent color stimuli as vectors in a three- 
dimensional LMS space, and the simulation algorithm is expressed in terms of transformations of this space. 
The algorithm replaces each stimulus by its projection onto a reduced stimulus surface. This surface is de- 
fined by a neutral axis and by the LMS locations of those monochromatic stimuli that are perceived as the 
same hue by normal trichromats and a given type of dichromat. These monochromatic stimuli were a yellow 
of 575 nm and a blue of 475 nm for the protan and deutan simulations, and a red of 660 nm and a blue-green 
of 485 nm for the tritan simulation. The operation of the algorithm is demonstrated with a mosaic of square 
color patches. A protanope and a deuteranope accepted the match between the original and the appropriate 
image, confirming that the reduction is colorimetrically accurate. Although we can never be certain of anoth- 
er's sensations, the simulation provides a means of quantifying and illustrating the residual color information 
available to dichromats in any digitized image. © 1997 Optical Society of America [S0740-3232(97)0 18 10-31 



1. INTRODUCTION 

Normal color vision is trichromatic. It is initiated by the 
absorption of photons in three classes of cones, the peak 
sensitivities of which lie in the long-wavelength (L), 
middle-wavelength (M), and short-wavelength (S) re- 
gions of the spectrum. 1 - 2 Therefore any color stimulus 
can be specified by three numbers, the cone responses; 
and all colors visible to the color-normal observer are in- 
cluded in a three-dimensional color space. Reduced 
forms of color vision arise from the effective absence of 
one of the retinal photopigments, 3 of the L type in protan- 
opes, the M type in deuteranopes, and the S type in 
tritanopes. 4 ' 5 For dichromatic observers any color stimu- 
lus initiates only two cone responses, and all colors that 
they can discriminate are included in a two-dimensional 
color space. Compared with trichromatic vision, dichro- 
matic vision entails a loss of discrimination and results in 
a reduced color gamut. In this paper we present a com- 
puterized method that allows color-normal observers to 
appreciate the range of colors experienced by dichromats. 
Designers of visual displays could use this method to 
simulate how their work would be seen by color-blind ob- 
servers. 

Early attempts to simulate the sensations of color-blind 
observers 6 " 8 faced one or other of two difficulties. On the 
one hand, one could apply an analog transformation to 



the whole image by using filters or eliminating one of the 
primaries of a three-color reproduction system: Such 
methods could be automated but would give an inaccurate 
simulation, since the spectral overlap of the cone sensitiv- 
ity curves means that removing a primary or removing a 
part of the spectrum is not equivalent to removing one of 
the classes of cone signal. On the other hand,, one could 
compute a separate colorimetric transformation for each 
element in the picture: With older technologies such a 
method would be cumbersome to apply to images of any 
complexity. Today, with the tools of digital imaging, and 
one of the sets . of cone fundamentals that are now 
available, 9-11 it is possible to use an automatic algorithm 
to simulate for the normal observer the appearance to the 
dichromat of any digitized image. 

We propose a computerized simulation of dichromatic 
vision in which the dichromat's color confusions and color 
palette are correctly represented for normal trichromats. 
Because dichromatic vision is a reduction of trichromatic 
vision, it should be possible to simulate for trichromatic 
observers the color gamut of dichromats. The simula- 
tion, for the normal trichromat, of dichromatic vision 
should not only reproduce the dichromat's color confu- 
sions but also offer a plausible color appearance. 

The color confusions are easily deduced from colorim- 
etry. Dichromats confuse lights that differ only in the 
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excitation of the missing class of cones and can discrimi- 
nate colors only on the basis of the responses of the two 
remaining cone types. In other words, one of the three 
colorimetric components that are required to specify a 
color stimulus in terms of the L- y Af-, and S-cone spectral 
responses is physiologically undetermined for a dichro- 
mat. 

As the dichromat does not see any change in the physi- 
ologically undetermined component, its value can be de- 
liberately chosen so as to imitate for the normal observer 
the appearance of colors for the dichromat. To this end, 
we make use of the color percepts of dichromats as re- 
ported in the literature. On the basis of these reports, we 
make several assumptions concerning the hues that ap- 
pear the same to dichromatic and to normal observers. 
These common hues define the color stimuli that remain 
invariant through the color-space transformations of our 
simulation. First, we assume that neutrals for normals 
are perceived as neutrals for dichromats. 12 Accordingly, 
no neutral stimuli are to be changed by the simulation. 
Second, we infer from reports on unilateral inherited 
color vision deficiencies that a stimulus of 575 nm is per- 
- ceived as the same yellow, and a stimulus of 475 nm as 
the same blue, by trichromats as by protanopes and 
deuteranopes. 12,13 Third, drawing upon a case of unilat- 
eral acquired tritanopia, we assume that the correspond- 
ing two hues for a tritanope are a red with a dominant 
wavelength of 660 nm and a blue-green with a dominant 
wavelength of approximately 485 nm. 14 

Taken together, these studies on the color perception of 
dichromats define the algorithm by which we determine 
the replacement value of the physiologically ineffective 
LMS component. In a preliminary paper, we have illus- 
trated the application of this algorithm to a natural scene 
of flowers. 15 Here we describe our algorithm in detail, 
and we demonstrate the application of the algorithm to a 
mosaic of randomly chosen color patches. 



2. METHODS 

To make our algorithm explicit, we represent color 
stimuli as vectors in a three-dimensional LMS space in 
which the orthogonal axes L, Af, and S represent the 
quantum catch for each of the three classes of cones. 
Such a representation allows the simulation algorithm to 
be expressed entirely in terms of LMS color-space trans- 
formations. 

A. Cone Fundamentals and LMS Color Space 
The LMS vector components (L Q , Mq , S Q ) of any given 
color stimulus Q are obtained from the spectral power dis- 
tribution of the stimulus, <p Q M (in radiometric units), 
and three spectral weighting functions /(X), m(X), and 
s(\) by numerical integration over the wavelengths \ of 
the visible spectral range: 

L Q = k j ¥> <? (X)7(X)dX, 
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Sq = k J <p Q (X)s(X)dX. (1) 

The spectral weighting functions 7(X), rh(\), and S(X) 
correspond to the action spectra of the three types of cone 
pigments as measured at the entrance of the eye. We 
used the cone fundamentals proposed by Stockman 
et aL lx in order to obtain /(X), m(X), and s(X); the detaUs 
of this procedure are described in Appc lix A. The com- 
mon scaling factor k was chosen such cnat Lq + Mq= 1 
for the White stimulus of maximum luminance obtainable 
on the color video monitor that we used for image presen- 
tation. 

B. Monitor Calibration 

We used a 20-in. color CRT monitor (Hitachi Model 
CM2086A3SG) driven by a computer graphics worksta- 
tion (Silicon Graphics Iris Indigo with GR2-XS24 graphics 
board; 3 x 8-bit color specification) for image display. 

We used a Bentham M300 spectroradiometer fitted 
with a 1200-line/mm diffraction grating (resolution: 2 
nm; wavelength step: 2 nm) and a secondary source of 
known spectral irradiance (traceable to the National 
Physical Laboratory, UK), to measure the spectral power 
distributions <p R M, <PgM> ***** of the red, green, 

and blue CRT primaries at the maximum intensity of 
each video channel. Table 1 lists the L» , the Mi , and the 
Si of the CRT primaries (i = R t G, B), which we ob- 
tained according to Eqs. (1) from the measured spectral 
power distributions. 

We used a Pritchard model 1980A-PL luminance meter 
to measure CRT luminances as a function of the digital 
video signal values and used these measurements to con- 
figure the computer graphics system for a linear relation 
between the digital pixel values and the CRT luminance 
for each of the R t G, and B video channels. 

C. LMS Color-Space Transformations 

According to the representation of color stimuli as tris- 
timulus vectors, additive color mixing on a color monitor 
can be written as 

Q = R Q V R + GqP g + BqP b , (2) 

where the vectors P R , P c , and P B represent the three 
CRT primaries at maximum intensity. The coefficients 
R Q , Gq> and B Q are weighting factors (ranging from 0 to 
1) that determine the relative contribution of each pri- 

Table 1, LMS Tristimulus Values for the Red, 
Green, and Blue Primaries and Nominal White for 
the Hitachi Model CM2086A3SG Monitor with the 
Respective Pixel Values at Maximum 0 





i = R 


i = G 


I = B 


i = W 


Li 


0.1992 


0.4112 


0.0742 


0.6846 


Mi 


0.0353 


0.2226 


0.0574 


0.3153 


Si 


0.0185 


0.1231 


1.3550 


1.4966 



°The L it M £ , and Si values were obtained according to Eqs. (1) from 
spectroradiometric measurements. 
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mary to the generated color stimulus Q. In terms of 
their LMS tristimulus values, Q and the three are ex- 
pressed as 



475 nm 



(i = R, G y B). (3) 



With the use of this representation, Eq. (2) becomes in 
matrix notation 



(4) 



The elements of the matrix T (L it "M it and S { for i 
= R, Q, B) are the LMS tristimulus values of the CRT 
primaries as given in Table 1. The inverse of Eq. (4), 
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V = T _1 Q, 



(5) 



allows-one to compute the pixel values Rq, Gq, and B Q 
for a color stimulus that is specified by L Q) M Q , and Sq . 
The practical application of our algorithm to digitized im- 
ages required three consecutive transformation steps for 
each pixel: (1) to compute the LMS specification Q from 
the original pixel values V by means of Eq. (4); (2) to ap- 
ply the simulation algorithm (Q-»Q'), which is described 
in the following subsection; and (3) to compute the result- 
ing pixel values according to Eq. (5), V = T Q'. 

D. Geometric Representation of the Algorithm 

Figure 1 shows a geometric representation of the algo- 
rithm. The orthogonal axes L, M, and S represent the 
light action on each of the three classes of cones, which 
have their maximal absorbance in the long-, middle-, and 
short-wavelength range of the visible spectrum, respec- 
tively. 

The limits of the parallelepiped define the gamut of 
color stimuli that could be obtained on the video monitor 
used for the simulation. W represents the nominal 
White stimulus obtained with Red, Green, and Blue video 
channel signals at maximum. E represents the brightest 
possible metamer of an equal-energy stimulus on this 
monitor. We have selected a metamer of the equal- 
energy stimulus E as the neutral color. OE represents 
the neutral stimuli for dichromats as well as for normal 
trichromatic observers. 

All stimuli that are perceived as neutral by dichromatic 
and normal trichromatic neutrally adapted observers lie 
between the origin O of the LMS space and the stimulus 

E. The neutral axis OE divides the surface of reduced 
stimuli into two half-planes, each of which is anchored on 
a point specifying an invariant hue for a given type of 
dichromat. The intersection of the surface of reduced 
stimuli with the parallelepiped representing the colors 
achievable on a CRT consists of two wings characteristic 
of protanopia and deuteranopia or two wings characteris- 
tic of tritanopia. The wings from OE toward the 475-nm 
and 575-nm locations represent the reduced stimuli sur- 
face for protanopic and deuteranopic stimulation [Fig. 




575 nm 
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Fig. 1. Geometric representation of the algorithm to. simulate 
dichromatic vision. The axes L, M, and S represent the quan- 
tum catch for each of the three classes of cones. The limits of 
the parallelepiped define the gamut of color stimuli that could be 
obtained on the video monitor used for the simulation. W rep- 
resents the nominal White stimulus obtained with Red, Green, 
and Blue video channel signals at maximum. E represents the 
brightest possible metamer of an equal-energy stimulus on this 
monitor. The line OE represents the neutral stimuli for dichro- 
mats as well as for normal trichromatic observers, (a) The 
wings from OE toward the 475-nm and 575-nm locations repre- 
sent the reduced stimuli surface for protanopic and deuteranopic 
simulation; (b) the wings from OE toward 485 nm and 660 nm 
represent the reduced stimuli surface for tritanopic simulation. 
The simulation algorithm consists in replacing the physiologi- 
cally undetermined component by the corresponding values on 
the respective reduced stimuli surface. This is illustrated for a 
given color stimulus Q, which is replaced by for the protan- 
ope, by % for the deuteranope, and by Q; for the tritanope. 
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1(a)], and the wings from OE toward 485 nm and 660 nm 
represent the reduced stimuli surface for tritanopic simu- 
lation [Fig. Kb)]. Given a stimulus Q by its LMS speci- 
fication, the algorithm replaces the undetermined compo- 
nent for a dichromat by the value corresponding to the 
projection of Q onto the wing, parallel to the direction of 
the missing fundamental axis. 

With the use of vector algebra, the equation for any 
stimulus Q' on a plane defined by the stimuli E, the 
monochromatic anchor stimulus A, and the origin O 
reads as 

(E X A)Q' = 0, (6) 

since the vector Q' is always orthogonal to the normal 
vector (E X A) of the plane. For given E 
= (L By M Ef S E ) and A = (L A , M A , S A ), this is a lin- 
ear equation of the coordinates Lq* , Mq> , and Sq* of 
stimulus Q': 

clLq* + bM Q , + cS Q , = 0, (7) 

with 

a = M E S A - S E M A , 

b = S E L A - L E S A , 

c = L E M A - M E L A . (8) 

Equation (7) with Eqs. (8) allows one to determine the 
replacement value of the physiologically ineffective LMS 
component. Thus the simulation algorithm (Q »-* Q') 
explicitly reads as, for the protanopic simulation, 

L Qr = ~{bM Q + cS Q )/a, 
M Q . = M Q , 

Sq* = Sq , (9) 

for the deuteranopic simulation, 
Lq* = Lq , 

M Q . = —(clLq + cS Q )/b, 
S Q , = S Q ; (10) 
and for the tritanopic simulation, 
Lqi = Lq , 

Mq* = Mq , 

S Q . = -(aL Q + bM Q )/c. 

As the two wings intersect at an angle along the neu- 
tral axis OE, the choice of the wing onto which a given 
stimulus Q is projected depends on the position of Q in 
relation to the neutral axis OE. The rules for this pro- 
jection are shown below: 



where the wavelength \ A indicates the monochromatic 
anchor stimulus of the wing onto which Q is projected. 

Figure 2 shows how the simulation is limited by the 
gamut of colors on the video monitor. All CRT colors are 
included in a parallelepiped in the LMS colorimetric 
space. If we project, parallel to one fundamental axis, 
the total volume of CRT colors onto a half-plane including 
the achromatic axis OE, we draw a polygon, and a part of 
the projection lies exterior to the shaded wings. The 
three parts of Fig. 2 show this type of parallel projection 
for each of the three forms of dichromatic vision: projec- 
tion along the L axis as for the protanopic simulation [Fig. 
2(a)], along the M axis as for the deuteranopic simulation 
[Fig. 2(b)], and along the S axis as for the tritanopic simu- 
lation [Fig. 2(c)], We have overcome the gamut limita- 
tion by starting the transformation with an image con- 
sisting of a subset of the colors obtainable in the RGB 
space of the monitor. . 

The original image data for the mosaic of square color 
patches [Fig. 3(a)] were obtained computationally by 
means of a pseudo-random-number generator under the 
constraint that all three projections onto the reduced 
stimuli surfaces were inside the gamut of the video moni- . 
tor (represented as shaded wings in Figs. 1 and 2). 

The original image and one of the simulated versions, 
each of which subtended 8 deg x 8 deg, were presented 
simultaneously on a 40-deg X 32-deg uniform neutral 
background (37 cd/m 2 , chromaticity of E). Dichromatic 
observers viewed the images (30- cd/m 2 average lumi- ; 
nance) at a distance of 0.5 m, and their task was to com-1 
ment on any noticeable difference. 

E. Observers 

Color-deficient subjects were classified on the basis of 
their performances on a battery of standard tests that in- 
cluded the Ishihara plates, the Panel D15, the desatu- 
rated D15, the Farnsworth-Munsell 100-hue, and 
matches with the Nagel anomaloscope. Two dichromatic 
subjects, GA and KBS, were selected to participate in all 
verification procedures. GA and KBS accepted matches 
in the whole red-green range with the Nagel anomalo- 
scope. For the protanope GA, settings of the brightness 
of the yellow field ranged from 30 when the Nagel red- 
green mixture was set at 0 to 3.5 when the mixture field 
was set at 73; for the deuteranope KBS, the correspond- 
ing values were 14 and 16.5 (a red -green mixture setting 
of 0 represents pure green, whereas 73 corresponds to 
pure red). A tritanopic subject was not available. 



3. RESULTS 

Figure 3 presets our simulation of the reduced color 
gamut seen by each class of dichromat, after transforma- 
tion of each pixel of the original mosaic image. The simu- 
lation of protanopic [Fig. 3(b)] and deuteranopic [Fig. 3(c)] 
vision is presented with the same hues: a blue and a yel- 



protanopic simulation: if 
deuteranopic simulation: if 
tritanopic simulation: if 



S Q /M Q < S e /M e , then k A 
Sq/Lq < S E /L Et then \ A 
Mq/Lq < M E IL E , then \ A 



575 nm; else k A - 475 nm; 
575 nm; else \ A - 475 nm; 
660 nm; else \ A = 485 nm, 
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(b) 



low. Differences in lightness are clearly visible. The 
small elements that appear red to the normal observer 
are dark for the protanope and light for the deuteranope. 
Differences in saturation are also visible. More satu- 
rated blue elements are present in the protanope simu- 
lation, and more saturated yellow elements are present in 
the deuteranope simulation. 

Figure 3 also shows that two elements that appear pur- 
plish red and orange-red to a normal observer may ap- 
pear, respectively, bluish and yellowish to a deuteranope. 
Consider, for example, the pink and orange-red elements 
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Fig. 2. Parallel projections of the video monitor gamut limits 
and reduced stimuli surfaces along (a) the L axis as for the pro- 
tanopic simulation, (b) the M axis as for the deuteranopic simu- 
lation, and (c) the S axis as for the tritanopic simulation. The 
shaded wings in (a) and (b) go toward 475 nm and 575 nm, and 
the wings in (c) go toward 485 nm and 660 nm. 



in line 6, columns 5 and 6 of the normal mosaic; their 
counterparts, in the deuteranopic mosaic are bluish and 
yellowish and more obviously different [Figs. 3(a) and 
3(c)]. The pink element is close in color to the flower of 
Pelargonium zonale, which the deuteranope John Dal ton 
judged "sky-blue" by daylight. 16 - 17 Dalton himself re- 
marks: "Red and scarlet form a genus with me totally 
different from pink. My idea of red I obtain from vermil- 
ion, minium, sealing wax, wafers, a soldier's uniform, 
etc.* 

The protanope GA and the deuteranope KBS each ac- 
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cepted the match between the original image and the 
transformation corresponding to his respective deficiency 
and could discriminate between the original and the other 
transformations. Careful inspection of every square of 
the mosaics paired for his deficiency did not help each ob- 
server in discriminating between the original and the 
simulated image. 



4. DISCUSSION 

A. Advantage of LMS Colorimetry 

Contrary to the CIE XYZ specification, which hides the 
relationships among the cone responses, the LMS specifi- 
cation, which decorrelates cone responses, is a very pow- 
erful colorimetric tool. In the case of the protan simula- 
tion, for example, it would not be enough to change the 
CIE X component, for any change in X affects both long- 
and middle-wavelength cones. Rather, we can derive a 
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secure algorithm from a three-dimensional vector space 
whose orthogonal ordinates correspond to the LMS 



B. Color Appearance 

As the simulation is based on color confusions exhibited 
by dichromats, the latter should not notice any difference 
between the original image and the simulated image cor- 
responding to their deficiency. Indeed, a protanope and a 
deuteranope were satisfied with the matches between the 
transformed image and the original image presented on 
the calibrated monitor screen used at the laboratory. 
The responses of dichromats are expected from classical 
measurements of dichromatic color matching 20 and vali- 
date our colorimetric stransformation. However, for any 
given image, there are an infinite number of alternative 
reductions that are metameric images for a particular 
type of dichromat. Only one of these reductions correctly 
demonstrates the color gamut enjoyed by dichromats. Tc 




(d> (c) 

Fig. 3. Reproduction of the video monitor display demonstrating the application of the simulation algorithm to a mosaic of randomly 
chosen color patches. Part (a) is the original, and simulations are given of how (a) is seen by (b) a protanope, (c) a deuteranope, and (d) 
a tritanope. (A colorimetrically exact reproduction of the video display cannot be guaranteed in the printed version.) 
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which extent does the final result demonstrate to the 
color-normal observer the color vision of the dichromats? 
We have followed unilateral dichromacy to select the 
wavelengths at which to anchor the reduced stimuli 
planes. Our choice lies in the range reported in analyses 
of hue matching or of color matching 12 ' 13 by unilateral 
dichromats. 

But we would make cautionary comments. First, the 
genetic basis for unilateral inherited dichromacy is not 
well understood. In the case of females, such observers 
could be*mterpreted as manifesting heterozygotes, whose 
normal X chromosome happened to be inactivated by ly- 
onisation in all the cones of one eye, 21 but such an expla- 
nation would not be applicable to males. In fact, there 
are few cases of unilateral dichromacy in which the good 
eye is convincingly normal: Often the Ishihara test is 
failed with the good eye. 22 %k 

A second and interesting problem in relying on unilat- 
eral dichromats is that the brain is plastic and seeks to 
minimize discrepancies between its inputs, as has been 
classically observed when the input to one eye is optically 
distorted, 23 such a process may minimize the discrepan- 
cies between the sensations evoked by the two eyes of the 
unilaterally color blind. 24,25 

A third cautionary comment concerns our implicit as- 
sumption that for a particular type of dichromat only one 
color subsystem is affected. Even though this assump- 
tion appears reasonable on the basis of the genetic expla- 
nation of dichromacy, there is also some evidence to the 
contrary. Chapanis 26 found the colors of red-green 
dichromats to be less saturated across the whole spec- 
trum, and more recently, Regan et al 21 reported that pro- 
tanopes and deuteranopes, as a population, had signifi- 
cantly higher tritan thresholds than normal trichromats. 



C. Abney Effect 

An ancillary advantage of our choice of hues to simulate 
the dichromats reduced sensation is that the Abney effect 
is minimized. The Abney effect, in which a variation of 
chromatic purity results in a hue shift, is absent for a par- 
ticular yellow and a particular bluish purple. 23 - 29 At 
575-nm dominant wavelength, as in our protan and deu- 
tan simulations, there would be hardly any hue change as 
the purity of the stimulus changes. At 475 nm and 485 
nm, the Abney effect, though present, remains small, 
within the gamut of chromaticities achievable by color 
mixture on a video screen. In fact, the unique hue locus 
is "curved, but mainly when it approaches the spectral 
locus. 30 Accordingly, stimuli within a reduced stimulus 
plane look uniform in hue to a color-normal observer, as 
we assume that they look to a dichromat. The choice of 
660 nm for the tritanope does not have the advantage of 
minimizing the Abney effect, but our first concern was to 
be consistent with a well-studied case of unilateral 
tritanopia. 14 

Actually, if we shift the wavelength at which the re- 
duced stimulus plane is anchored, from 575 nm to 570 nm 
for instance, the single hue sensation for the normal ob- 
server disappears. This hue breakdown is manifest with 
natural images in which many spatial frequencies are 
included. 31 A possible explanation of this phenomenon 
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might be that the size of the Abney effect depends not 
only on the dominant wavelength but also on spatial 
frequency. 32 

D. Contrast 

Finally, beyond the geometric justification, it is worth 
wondering whether the fact that the projection algorithm 
required the limitation of the original image colors to a 
subset of CRT colors implies that a dichromat will see, in 
the case of a small range of chromaticities, colors that cor- 
respond to a more extended color gamut than that for nor- 
mal observers. Such a phenomenon is unexpected. 

E. Polymorphism of Dichromatic and Normal Vision 
Since dichromats of a given type are known to vary 
slightly in their spectral sensitivity, we might expect 
them to vary in their willingness to accept a match be- * 
tween the original image and the appropriate reduction. 
Moreover, polymorphism of visual pigments has been 
demonstrated within the color-normal population as well 
as within classes of dichromats, 33 " 38 and the available 
fundamentals represent only an average observer. In a 
recent study - using the electroretinogram, Neitz et a/. 39 
have found four distinctly different spectral sensitivities 
of the L- and M-cone pigments among 12 dichromatic ob-; 
servers. The M-cone pigment in seven of eight pro tan- 
opes had a spectral peak of 530 nm, and the pigment of 
the eighth had a peak at 537 nm. The peak sensitivity of 
the L-cone pigment was at 558 nm in two of four deutera- 
nopes and at 563 nm in the other two. The 5-7-nm shifts 
within each class of cone pigments were accounted for by- 
having alanine or serine at site .180 of the amino-acid se- 
quence of the pigment molecule. 39 The serine/alanine 
polymorphism has been shown earlier to be responsible 
for variations in color matching among color-normal 
observers. 34 Of course, refinements of our algorithm 
could simulate alternative forms of, say, deuteranopia for 
alternative types of normal observers. 



5. CONCLUSION 

The availability of computer-controlled color displays has 
allowed us to develop an algorithm to simulate for normal 
observers how the dichromat perceives a complex colored 
scene. Although one can never know for certain the qual- 
ity of another's sensations, our simulation allows us to 
quantify the range of residual color information available 
to dichromats in any digitized image. 



APPENDIX A: CONE FUNDAMENTALS 

The computational procedure to obtain the spectral cone 
contribution functions m(\), and s(k) for Eqs. (1) 

was as follows: We used the values of columns 5-7 of. 
Table 8 in Stockman et al. 11 as c5(X), c6(\), and c7(X): 

/'(A) * 0.68273 x 10 c5a >, 

m'(\) = 0.35235 x Xf<* K \ 

s'(X) = 1.00000 x lO 07 ^. (Al) 
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These are 2-deg cone fundamentals derived from the 
color-matching functions of the CIE 1964 supplementary 
standard colorimetric observer. The weighting factors 
0.68273 and 0.35235 in Eqs. (Al) are those indicated in 
Ref. 11 in order to obtain an approximation to the CIE 
modified V x function V M (\) by adding I'M to m'(\). 
The use of these weighting factors allows for convenient 
scaling: 

7(\) = Z'OO J 2 P'OO + *»'<ML 

m(\) = m'(\) / 2 [I'M + m'(X)], 

i(X)-V(X) / 2 s'(X), r (A2) 

with the result that the nominal equal-energy stimulus 
then has L + M = 1 and S = 1. 
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